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Abstract 

We study the effects of dark-matter annihilation during the epoch of big-bang nucle¬ 
osynthesis on the primordial abundances of light elements. We improve the calculation 
of the light-element abundances by taking into account the effects of anti-nucleons 
emitted by the annihilation of dark matter and the interconversion reactions of neu¬ 
tron and proton at inelastic scatterings of energetic nucleons. Comparing the theoreti¬ 
cal prediction of the primordial light-element abundances with the latest observational 
constraints, we derive upper bounds on the dark-matter pair-annihilation cross section. 
Implication to some of particle-physics models are also discussed. 


Recent cosmological and astropliysical observations have revealed that about 26% of the 
mass density of the present universe is occupied by dark matter [1]. This fact suggests that 
there exists a stable (or very long-lived) particle or held which behaves as a non-relativistic 
object in the present universe. However, the particle-physics nature of dark matter is almost 
unknown yet, and it is an important task to acquire information about it. 

From particle-physics point of view, a stable particle has attracted many attentions as a 
candidate for dark matter; in our analysis, we assume that some stable particle plays the role 
of dark matter and denote it as X. In order for an efficient production of the dark-matter 
particle in the early universe, it is often the case that X has interactions with some of the 
standard-model particles. In such a case, we may derive constraints on the properties of 
dark matter by studying the effects of dark-matter pair annihilation into standard-model 
particles in the early universe. 

One of the important effects of dark-matter pair-annihilation is on the big-bang nucle¬ 
osynthesis (BBN). If dark matter pair-annihilates into charged particles during or after the 
BBN, they induce photodissociation processes of light elements synthesized via the BBN 
reactions. In addition, energetic hadrons are often produced as a consequence of the pair- 
annihilation; if so, hadrodissociation of the light elements are induced. Because the standard 
BBN scenario predicts the light-element abundances which are more-or-less consistent with 
observations, too large pair annihilation cross section is excluded. Indeed, such a constraint 
has been intensively discussed in literatures [aslIllElEllTllH]® 

Recently, there have been progresses in the observational determination of the primordial 
abundances of the light elements. In particular, uncertainty in the primordial deuterium (D) 
abundance has been signihcantly reduced. Such a progress has a large impact on the BBN 
bounds on the dark-matter properties. 

In this letter, we revisit the BBN constraint on the annihilation cross section of dark 
matter, taking into account the recent progresses in the observation of the light-element 
abundances. For a reliable calculation of the light-element abundances, we improve the 
treatment of hadrodissociation; in particular, we have newly included the effects of anti¬ 
nucleons emitted by the annihilation of dark matter as well as the effects of interconversion 
between neutron and proton at inelastic scatterings of energetic nucleons. Then, comparing 
the theoretical predictions with observational constraints on the light-element abundances, 
we derive upper bounds on the annihilation cross section of dark matter. 

We hrst summarize the current status of the observational constraints on the primordial 
abundances of light elements. The primordial abundance of D is inferred from D absorption 
in damped Lya systems (DBAs). Recently Cook et al. [lU] observed a DBA toward QSO 
SDSS J1358-I-6522 and performed very precise measurement of D. They also reanalyzed other 
four previously known DBAs and, using the total hve DBA samples, obtained the primordial 
D abundance as 

(D/H)p = (2.53 ± 0.04) x 10-^ (1) 

"^^Besides photodissociation and hadrodissociation, annihilation of light dark matter {mx ^ 0(10) MeV) 
affects the BBN by changing the temperature ratio between neutrinos and photons, from which constraints 
are obtained [9]. 
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where (A/B) denotes the ratio of number densities of light elements A and B, and p indicates 
the primordial value. Notice that the error is smaller by a factor of 5 than that adopted 
in our previous study [6]. This progress in the measurement of D leads to more stringent 
constraints on dark-matter annihilation as seen later. 

As for the primordial mass fraction, Yp, of helium 4 (^He), a new determination with the 
use of the infrared as well as visible ^He emission lines in 45 extragalactic HII regions was 
reported in Ref. where Yp = 0.2551 ± 0.0022 is obtained. More recently, Aver, Olive 
and Skillman [12] reanalyzed the data of Ref. m and estimated the ^He abundance using 
Markov chain Monte-Carlo analysis. They obtained 


Yp = 0.2449 ± 0.0040, 


( 2 ) 


which we adopt in this letter. 

We also use a constraint on ^He/D which is derived from D and ^He abundances observed 
in protosolar clouds [ 13 ]; taking into account that the ratio ^He/D increases monotonically 
in time, we adopt 

(^He/D)p < 0.83 + 0.27. (3) 

This observational constraint is the same adopted in Ref. [6]. 

In the previous studies, constraint based on the lithium 7 (^Li) abundance was also 
discussed. However, the situation of the ^Li observation is now confusing. The observed ^Li 
abundances in metal-poor halo stars showed almost a constant value (log^g(^Li/H) ~ —9.8) 
called Spite plateau which was considered as primordial. However, the recent observation 
found much smaller ^Li abundances (logio(^Li/H) < —10) for more metal-poor stars [T4] . 
Since we do not know any mechanism to explain such small abundances, we do not use ^Li 
to constrain the properties of dark matter in this letter. We do not use ®Li either because 
®Li abundance is observed as the ratio to the number density of ^Li. 

In order to derive constraints on the dark-matter properties, we calculate the primordial 
abundances of the light elements, taking into account the effects of dark-matter annihilation. 
Our calculation of the light-element abundances is based on Refs. [ISile] with the modi- 
hcations explained below. The Boltzmann equations for the evolution of the light-element 
abundances have the following form: 


dUA 

dt 


SHha = 


dUA 

+ 

1 

■ 

-f 

dnA 

+ 

1- 

■ 

dt 

SBBN 

dt 

photodis 

dt 

hadrodis 

dt 


( 4 ) 


J p-^n 


where ua denotes the number density of the light element A, and H is the expansion rate 
of the universe. Here, [duA/dt]sBBN denotes the effects of the standard BBN reactions while 
the other terms in the right-hand side are the effects of dark-matter annihilation, i.e., those 
of photodissociation, hadrodissociation, and p ■(-)■ n conversion. The reaction rates due to 
the dark-matter annihilation are proportional to the annihilation rate of dark matter which 


is given b\l^^ 


rannihilation — , 


( 5 ) 


In a series of papers [muB], the effects of decaying particles had been studied. The Boltzmann equations 
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where nx is the number density of dark matter. In addition, {av) is the annihilation cross 
section, with which the Boltzmann equation for the evolution of the number density of dark 
matter is given by 


+ ^Hnx = (av) - n^) , (6) 

where nx,eq is the equilibrium value of the number density of dark matter. We assume that 
the dark-matter annihilation occurs through an s-wave process so that (av) is independent 
of the relative velocity of dark-matter particles in the non-relativistic limit. The effects of 
the pair annihilation do not change rix signihcantly during and after the BBN epoch because 
those epochs are long after the freeze-out time of dark matter. Then we use the following 
time-evolution of nx'- 


mx V «o / ’ 

where Mpi ~ 2.4 x 10^® GeV is the reduced Planck scale, Hq is the Hubble constant, 
and mx are the density parameter and the mass of dark matter, respectively, and a{t) and 
Oq are the scale factor at the cosmic time t and at present, respectively. (In our numerical 
calculation, we use Hq = 68 km/sec/Mpc and flxh‘^ = 0.12.) [dn^/dt]photodis, [<^ri^/<^t]hadrodis, 
and [dnA/dt]p^ri are proportional to n\{av), and the effects of the dark-matter annihilation 
become more efficient as the annihilation cross section {av) increases. In order not to affect 
the light-element abundances too much, {av) is bounded from above. 

Next, we summarize the new points in the calculation of the light element abundances 
compared to Refs. [lade]: revision of the SBBN reaction rates and the treatment of the 
hadrodissociations. In the present study, we renewed reaction rates by adopting the results 
of Ref. [I?] @ In order to take into account the uncertainties in the reaction rates, we use the 
Monte-Carlo simulation to estimate the errors of the light element abundances by assuming 
that the errors of the reaction rates obey Gaussian distributions. 

As for the hadrodissociations, we have improved the treatment of the hadronic showers 
initiated by injections of energetic hadrons into the thermal plasma, including the following 
effects!^ 

1. Effects of anti-nucleons emitted from dark-matter annihilation. 

used in the present analysis can be obtained from those in Ref. [151 [1^ by replacing Fdecay Fannihiiation, 
where Fdecay is the decay rate defined in Refs. [HI [16], with properly rescaling the numbers of final-state 
particles produced by the pair annihilation of dark matter using the fact that two dark matters participate 
in the pair-annihilation process (instead of one for the decay process). 

For another recent study of nuclear reaction rates, see also Ref. [18], which we do not use in our analysis 
because there is a technical difficulty in implementing the errors given in Ref. [18] into our Monte-Carlo 
analysis. 

"^■^For more details, see Ref. [T9] . 


nx{t) 


( 7 ) 
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2. Effects of interconversion reactions between (anti-) neutron and (anti-) proton, with 
which injected and secondary-produced beam nucleons, as well as target nucleons, 
change their charges at the time of the inelastic scattering. 

Concerning the effects of anti-nucleons, we have considered scatterings of anti-nucleons off the 
background protons and ^He’s. The scatterings of anti-nucleons off the background protons 
produce high-energy nucleons, which can destroy ^He’s and produce copious high-energy 
daughter particles (i.e., nucleons and light elements). Although the inelastic scatterings of 
anti-nucleons off ^He’s produce high-energy daughter particles as well, we have neglected all 
the hnal-state particles of the scatterings in the subsequent calculation because of the lack of 
sufficient experimental data for those inelastic scatterings. It is expected that the inclusion 
of such hnal-state particles would only make the BBN constraints severer, and hence our 
treatment gives conservative bounds. 

With the inclusion of the effects of anti-nucleons, the constraints become stronger by 10% 
- 30% for mx ~10 GeV - 1 TeV. On the other hand, by the latter effects (i.e., the effects 
of the interconversions at inelastic scatterings), the constraints become weaker by 50% - 
80%. This is due to the fact that a high-energy proton, which is interconverted from a beam 
neutron, tends to be stopped more easily through electromagnetic interactions in the plasma 
than the neutron with the same energy. 

Now, we are at the position to discuss the bounds on the annihilation cross section of 
dark matter (crn). The effects of the dark-matter annihilation on the BBN depend on the 
hnal-state particles produced by the annihilation process. Here, we derive the upper bound 
on (an), assuming that the annihilation process is dominated by one of the following modes: 

• XA ^ W+W-, 

• XX —)■ qq (with q = u,d, s, c, b). 

We consider only the cases where the pair annihilation results in the production of two par¬ 
ticles with an identical mass. Thus, the energy of individual hnal-state particles is equal to 
mx- For the calculation of the photodissociation rate, it is necessary to acquire the total 
amount of the energy injection in the form of electromagnetic particles due to the pair anni¬ 
hilation. In addition, the study of the hadrodissociation and the p n conversion processes 
requires energy distributions of hadrons (in particular, proton, neutron, and pions) produced 
by the annihilation of dark matter. The decay, cascade, and hadronization processes of the 
standard-model particles are studied by using PYTHIA 8.2 package [2U] . 

Comparing the theoretical predictions of light-element abundances with observational 
constraints, we derive the bounds on the annihilation cross section (crv). The results are 
shown in Figs. [Hand [2] for W~^W~ and qq emissions, respectively!^ On the figures, only 
the constraints from D and ^He are shown because that from ^He/D is too weak to show up 
in the figures. We can see that the constraints from D are much more stringent than those 
from ^He. 

have also performed an analysis for the case where the dark matter dominantly annihilates into ZZ. 
We have found that the bound is almost the same as the case with the W~^W~ final state. 
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Figure 1: Upper bound at 95% C.L. on the annihilation cross section into the W~^W~ mode 
as a function of the dark-matter mass. The upper bound comes from the observational 
deuterium abundance. We also plot the theoretical prediction of the annihilation cross 
section for a pair of Winos into the W~^W~ mode taken from Ref. [5]. 

The upper bounds on (av) can be converted to constraints on model parameters once 
the particle-physics model of dark matter is specified. First, let us consider the case where 
the neutral Wino is dark matter in a supersymmetric model. The neutral Wino is the 
superpartner of the neutral SU{2)l gauge boson and one of the well-motivated candidates of 
dark matter, which dominantly pair-annihilates into a W^W~ pair. In order for the thermal 
relic Wino to be dark matter, the Wino mass is required to be ~ 3 TeV [21]. However, Wino 
with smaller mass can also account for dark matter if it is non-thermally produced [2^ |23] . 
In such a case, a larger value of {av) is possible. In Fig. [H we also plot the cross section for 
the process W^W~ as a function of the dark-matter mass (i.e., the Wino mass). 

We can see that there exist mass ranges in which the cross section is larger than the upper 
bound obtained from the BBN. Assuming that neutral Wino is dark matter, the Wino mass 
is constrained to be 


320 GeV <M^< 2.3 TeV or > 2.5 TeV. (8) 

With the present analysis, the BBN bound on the Wino mass has become more stringent 
than that obtained in the previous BBN analysis. Neglecting the constraints from ®Li and 
^Li, the bound was given by 250 GeV < < 2.3 TeV or > 2.4 TeV [5]. We can 
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Figure 2: Upper bounds at 95% C.L. on the annihilation cross section into the qq mode (hn, 
dd, ss, cc, and bb) as a function of the dark-matter mass. The upper bounds come from the 
observational deuterium abundance. 

see that, in the Wino dark-matter scenario, the lowest possible value of is signihcantly 
increased. In addition, the scenario of the thermal Wino dark matter (with ~ 3 TeV) is 
allowed. Next, let us discuss the case of thermal relic dark matter, assuming that the dark 
matter dominantly annihilates into quark pairs. In such a case, (av) ~ 3 x 10“^® cm^/sec 
is required, which results in mx > 25 — 35 GeV. There are some remarks related with 
the recent reports of the possible y-ray line excess from the Galactic center at around GeV 
energies. In Ref. m, for example, it is claimed that the signal can be htted by the qq 
emission from annihilation of dark matter with the mass of 36 — 51 GeV@ Such a region is 
still consistent with the BBN bounds. 

Finally, let us comment on other collider and astrophysical constraints. One constraint 
is from the direct search of charged Wino at the LHG experiment. Because the charged and 
neutral Winos are almost mass-degenerate, the charged Wino becomes relatively long-lived 
(i.e., ~ 0(0.1 ns)) if the neutral Wino is the LSP [2Sl EZ]- Then, the charged Wino 

production may leave disappearing track at the LHG. The bound from the disappearing- 
track search is currently > 270 GeV at 95% G.L. [281129] . which is weaker than the 
BBN bound. 

Another constraint is estimated from the study of the energetic y-ray fluxes from Milky 
Way satellites. In some case, sizable y-ray fluxes are expected from those satellites, although 
no signal has been observed yet. The negative search for the y-ray signal results in a 
lower bound on the mass of dark matter. For example. Ref. [30] obtained the constraint of 

other efforts, see also Refs. [25]. 
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320 GeV < < 2.25 TeV or > 2.43 TeV for the Wino dark-matter scenario with the 

study of four nearby dwarf spheroidal galaxies (dSphs) whose dark-matter profiles can be 
obtained with stellar kinematic data. In addition, with extra assumptions to determine the 
dark-matter profiles, the Fermi-LAT collaboration EH obtained 1 TeV < <1.7 TeV or 

> 2.8 TeV. The Fermi-LAT collaboration also obtained mx ^ HO and 150 GeV for the 
thermal relic dark matter scenario with W^W~ and qq emissions, respectively. Although 
these constraints are comparable to or severer than the BBN ones, constraints from y-ray in 
general may suffer from uncertainties of the density prohles of dark matter in galaxies. 

The recent observation of the anti-proton flux in high-energy cosmic ray [32] also puts 
a bound on the mass of dark matter [33]. Importantly, however, the theoretical prediction 
of the anti-proton flux from the dark-matter annihilation has large uncertainty because the 
anti-proton flux is sensitive to the model of cosmic-ray propagation as well as the dark-matter 
density profile of the Milky Way. Using the result of Ref. [3l], for example, the anti-proton 
constraint on the Wino dark matter scenario is > 240 — 250 GeV. In addition, for the 
case of qq emission for the annihilation process, the mass of dark matter is constrained to 
be mx ^ 20 GeV, assuming the annihilation cross section to realize the thermal relic dark 
matter scenario. These constraints are weaker than the BBN constraints obtained by our 
analysis. (However, they may become more stringent if one takes propagation models or 
density profiles other than that adopted in Ref. [34].) 

Since dark-matter annihilation affects the ionization history of the universe, the observa¬ 
tions of the cosmic microwave background (GMB) constrain the annihilation cross section. 
The recent Planck obtained (crn) < 4 x 10“^® cm^s“^/“^(mx/GeV) [1], where / is the frac¬ 
tion of the injected energy that goes into the background plasma. / is expected to be 
~ 0(0.1); for example / ~ 0.15 for the annihilation into and / ~ 0.2 for the an¬ 

nihilation into ii~ or bb [T]. If we take / ~ 0.2 as a typical value, the Planck result 
leads to the constraints, 400 GeV < < 2.1 TeV or > 2.6 TeV for the Wino dark 

matter, and mx ^ 15 GeV for the thermal relic dark matter with qq emission. (For more 
accurate bounds, precise calculation of / is necessary for each model.) Therefore, the GMB 
constraints are comparable to the BBN constraints. 

In summary, we have studied effects of dark-matter annihilation in the early universe 
on the abundances of the light elements synthesized during the BBN epoch. If the dark- 
matter annihilation results in the production of electromagnetic and hadronic particles, they 
affect the abundances of the light elements through photodissociation, hadrodissociation, 
and p o n conversion processes. We have calculated the abundances of D, ^He, and "^He, 
taking into account above processes. In the study of the effects of dark-matter annihilation, 
we have improved the treatment of the hadrodissociation processes. In particular, (i) we have 
included the effects of the anti-nucleon emitted by the annihilation process of dark matter, 
and (ii) we take account of the interconversion reactions between neutron and proton at 
inelastic scatterings. Then, comparing the theoretical prediction with the latest observational 
constraints on the primordial abundances of the light elements, we have derived the upper 
bounds on the pair annihilation cross section of dark matter. We found that the latest 
results on the measurements of the D abundance, which result in a precise determination of 
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the primordial abundance of D, have a strong impact. For the case of Wino dark matter, 
for example, the BBN constraint requires the Wino mass to be 320 GeV < < 2.3 TeV 

or > 2.5 TeV. 
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